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We present a comparative analysis of superconducting and charge-density-wave orders in the spin- 
fluctuation scenario for the cuprates. That spin-fluctuation exchange gives rise to d-wave supercon¬ 
ductivity is well known. Several groups recently argued that the same spin-mediated interaction may 
also account for charge-density-wave order with momenta (Q,0) or (0, Q), detected in underdoped 
cuprates. This has been questioned on the basis that charge-density-wave channel mixes fermions 
from both nested and anti-nested regions on the Fermi surface, and fermions in the anti-nested 
region do not have a natural tendency to form a bound state, even if the interaction is attractive. 

We show that anti-nesting is not an obstacle for charge order, but to see this one needs to go be¬ 
yond the conventional Eliashberg approximation. We show that in the prefect nesting/antinesting 
case, when the velocities of hot fermions are either parallel or antiparallel, the onset temperatures 
in superconducting and charge-density-wave channels are of comparable strength for any magnetic 
correlation length £. The superconducting T sc is larger than Tcdw, but only numerically. When 
the velocities of hot fermions are not strictly parallel/antiparallel, T c a w progressively decreases as £ 
decreases and vanishes at some critical £. 


I. INTRODUCTION 

The experimental discovery of static charge-density- 
wave (CDW) order in underdoped cuprates^ - — has re¬ 
ignited theoretical studies of the mechanism of CDW in¬ 
stability and its interplay with d-wave superconductivity. 
The CDW order accounts for a number of properties of 
the pseudogap phase in the cuprates and the understand¬ 
ing of the mechanism of CDW instability is an essential 
step towards the understanding of the pseudogap. 

Several groups^ - — analyzed charge order with diago¬ 
nal momentum ( Q , ±Q), which we refer to as CDW-diag, 
and compared it with with superconductivity within the 
spin-fluctuation (hot spot) scenario. It has been ar¬ 
gued that the hot spot model has an approximate SU(2) 
particle-hole symmetry, which makes the onset temper¬ 
atures for CDW-diag and superconducting order almost 
equal for large magnetic correlation length. It has been 
propose d 10 ’ 11 that the pseudogap may be due to the 
fact that over a wide range of T the system cannot de¬ 
termine between near-degenerate superconducting and 
CDW-diag orders. 

It turned out, however, that the CDW order in the 
cuprates has momenta Q = ( Q , 0) or (0, Q) along X or 
Y directions in the momentum s pac e rather than along 
the diagonals. We will follow Ref. Il3l and refer to this or¬ 
der as CDW-x. The experimental value of Q for CDW-x 
order in the cuprates is close to the distance between hot 
spotsi Fermions in the hot regions are the ones which 
mostly participate in magnetically-mediated interaction, 
and the closeness of the experimental Q to the distance 
between hot spots fueled speculations that the same mag¬ 
netic fluctuations which favor charge order with diagonal 
momenta ( Q,Q ) may also be responsible for CDW-x in¬ 
stabilit y 11 ’ 14- — . 


The magnetic scenario for CDW-x order is appeal¬ 
ing by two reasons. First, the CDW-x order between 
hot fermions develops together with pair-density-wave 
(PDW) orde r 11 ’ 12 ’ 18- — , and the combination of the two 
explains specific features of ARPES dat a 15 ’ 17 ’ 19 ’ 23 . Sec¬ 
ond, a Ginzburg-Landau (GL) analysis shows 12 i 15 i 16 i 24 i 25 
that CDW-x order breaks not only U( 1) translational 
symmetry, as is expected for any incommensurate charge 
order with a complex order parameter, but also C 4 lat¬ 
tice rotational symmetry and time-reversal symmetry. 
The breaking of C 4 is the consequence of the fact that 
CDW/PDW order develops in the form of stripes, and 
the breaking of time-reversal symmetry is the conse¬ 
quence of the fact that CDW-x order parameters between 
the pairs of hot fermions with center of mass momentum 
k and k develop with relative phase ±7r/2. These two 
order parameters transform into each other under time 
reversal, and by selecting one of these states, the system 
spontaneously breaks time-reversal symmetry. Both C 4 
symmetry breaking and time-reversal symmetry breaking 
have been observed in the experiments^^ - —. 

The GL analysis assumes that CDW-x order does de¬ 
velop at a finite T c d w , and that this temperature is com¬ 
parable to T sc for d-wave superconducting (SC) insta¬ 
bility, which also develops between hot fermions due 
to magnetically-mediated interaction— The mean-field 
value of T sc is expected to be larger than T cc j w sim¬ 
ply because superconductivity necessarily involves pairs 
of fermions with opposite directions of Fermi-velocities 
(nesting), while CDW-x instability involves pairs of 
fermions whose Fermi velocities are generally at some 
finite angle with respect to each other. Still, if mean- 
field values of T sc and T c d w are comparable, the effects 
beyond mean-field (e.g., the pre-emptive breaking of dis¬ 
crete symmetries for CDW-x order) may lift T c( j w above 
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T sc . This calls for a detailed comparative analysis of 
the onset temperatures of SC and CDW-x orders at the 
mean-field level, by which we mean ladder approxima¬ 
tions for the corresponding vertices (see Fig. 2). 

In a recent paper—, the two of us compared the onset 
temperatures of SC and CDW-x instabilities by analyz¬ 
ing the structure of the diagrammatic series for SC and 
CDW-x vertices in the quantum-critical regime, when the 
magnetic correlation length £ is infinite and the fermionic 
self-energy E(w) has a non-Fermi liquid form y/u and ex¬ 
ceeds the bare u> term at small frequencies^. We consid¬ 
ered the generic case of some finite angle between Fermi 
velocities of hot fermions separated by (Q,0) or (0,Q) 
and found that kernels of the “gap” equations are log- 
arithmical in both SC and CDW-x channels, and the 
prefactor for the logarithm for CDW-x channel is only 
numerically smaller than the one for SC channel. We cut 
the logarithms by temperature and found that T sc and 
Tcdw are of comparable strength at infinite £. At a finite 
£, the logarithm in T in the CDW channel is cut already 
at T = 0, and, as a result, T c dw decreases and vanishes at 
some finite £. The SC instability, on the other hand, is 
not critically affected by decreasing £ and, in the absence 
of impurity scattering, survives at any £. 

Although this analysis is plausible, the cutting of the 
logarithm by temperature is not a rigorously justified 
procedure in the quantum-critical regime because the log¬ 
arithm in a particular cross-section of the ladder series is 
cut by T (or £ _1 ) only if we set the frequencies of exter¬ 
nal fermions ±w m to zero. At a finite co m the logarithms 
in the SC and CDW-x channels are already cut by oj m 
even at T = 0 and £ = oo, and one has to go beyond the 
leading logarithmical approximation to rigorously ana¬ 
lyze the emergence of SC and CDW-x (see below). 

In this communication we present the alternative anal¬ 
ysis of SC and CDW-x instabilities, which does not rely 
on a comparison of logarithms in the perturbation theory. 
Specifically, we show that at infinite £ the full linearized 
equation for CDW-x order parameter, from which one ex¬ 
tracts the temperature of CDW-x instability, differs from 
the corresponding equation for the SC order parameter 
only by the strength of the effective coupling. We use 
recent non-perturbative result s 33 ' 34 for SC instability in 
the quantum critical regime which show, among other 
things, that the instability develops at any value of the 
coupling, and relate T sc and T cc j w . To see this, we fo¬ 
cus on the seemingly “worst case scenario” for CDW-x 
instability, when the Fermi surface is horizontal or ver¬ 
tical in hot regions and CDW order involves a half of 
fermions in the nested region and a half in the anti¬ 
nested region (see Fig. 1). We note in passing that an 
almost nested/antinested Fermi surface agrees well with 
ARPES data for B^S^CaC^Os+j; (the ratio of antipar¬ 
allel/parallel velocities for hot spots 1 and 2 in Fig. 1 is 
13.6, same for parallel/antiparallel velocities at hot spots 
3 and 4, see Ref. 0 ). For a homogeneous d-wave SC 
instability, nesting/anti-nesting is not an issue because 
the two fermions in the particle-particle channel neces¬ 


sarily have momenta k and — k and, hence, the same efc. 
For CDW-x instability, the kernel of the “gap” equation 
in the nested region is same as in the SC channel, i.e., 
the product of the two Green’s function G’s has (without 
self-energy) the same form l/(w 2 + e|) as the product of 
the two G's in the SC channel. On the other hand, the 
product of the two G’s for CDW-x in the anti-nested re¬ 
gion gives 1 — efc) 2 , and this combination by itself 

does not lead to Cooper logarithm after integration over 
efe and Lu m because it contains a double pole as a func¬ 
tion of either oj m or e*,. However, the true kernel in the 
anti-nested region contains the product of two G’s in the 
combination with the bosonic propagator, and the latter 
also depends on efc and ui m and contains poles (as a func¬ 
tion of efe) and branch cuts (as a function of u> m ) in both 
half-planes of complex ek and u) m . It then becomes an is¬ 
sue whether the contribution from the poles/branch cuts 
in the bosonic propagator yields the result comparable in 
magnitude to the one in the nested region. 

To analyze this issue, we do the same trick as was re¬ 
cently used in the analysis of the optical conductivity 
in the cuprates^ and re-express the set of two coupled 
equations for CDW-x order parameters in the nested and 
anti-nested regions as the single equation for the CDW-x 
order parameter in the nested region with the effective 
interaction from a second-order composite process involv¬ 
ing fermions in the anti-nested region (see Fig. 2). The 
effective composite interaction y com involves two G’s [the 
ones whose product gives l/(iu> — e*,) 2 ] and two spin- 
fluctuation propagators x{^,lo). We evaluate the prod¬ 
uct by integrating over k and w m and compare Xcom with 
the interaction in the SC channel, which is a single spin- 
fluctuation propagator. We show that the effective in¬ 
teraction is comparable to the original y, both in the 
Fermi liquid regime at moderate £ and in the quantum- 
critical regime at large £. That the two interactions 
are comparable by magnitude may seem strange because 
Xcom would vanish if we approximated spin-fluctuation 
propagators by their values between the particles on the 
Fermi surface, as it is done in the Eliashberg approxi¬ 
mation. However, this approximation is rigorously jus¬ 
tified for electron-phonon interaction, for which correc¬ 
tions to Eliashberg approximation are small in the ratio 
of phonon velocity to Fermi velocity, while for electronic 
pairing mechanism it is justified only in the artificial limit 
of large number of fermionic flavors N (Refs. i, 0, 0- 
[38b . which we do not impose here. For the physical case 
of N = 1, there is no parameter which would allow one 
to neglect the dependence on e*, in the spin-fluctuation 
propagators. Keeping these dependencies, we find that 
the integral which determines Xcom is non-zero due to the 
poles in the two bosonic propagators, considered as func¬ 
tions of efe, and, moreover, its magnitude is comparable 
to the original spin-fluctuation propagator y- 

We analyze both Fermi-liquid and quantum-critical 
regimes and show explicitly that T c <j w for CDW-x order 
is non-zero and differs from SC T sc only by a numeri¬ 
cal factor which increases as magnetic correlation length 
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FIG. 1. The Brillouin zone, the magnetic Brillouin zone, and 
FS using the same dispersion as in Ref. [h|. hot spots 1,2 
and 3,4 are defined as points on the Fermi surface that inter¬ 
sects with the magnetic Brillouin zone. The red wavy lines 
represents the interaction mediated by spin fluctuations of 
momentum (7r,7r). Hot spots -1, -2, -3, -4 (not shown) have 
momenta opposite to those of hot spots 1,2,3,4, respectively. 


£ gets larger. We then analyze the effect of deviation 
from perfect nesting/antinesting and show that T cc j w gets 
progressively reduced upon decreasing £ and eventually 
vanishes above some critical correlation length. 


II. THE MODEL 

We consider the same model as in earlier stud- 
ie 3 & i 10 i 32 i 36 i 37 i 39 - ~— : fermions in hot regions, interact¬ 
ing by exchanging Landau-overdamped magnetic fluctu¬ 
ations peaked at K = ( 7 r, 7 r). For the bulk of the paper 
we assume that the Fermi velocities of hot fermions sep¬ 
arated by {Q, 0) or (0, Q) are either parallel or antipar¬ 
allel. We approximate the dispersion of hot fermions by 

= Vyky and efc = —v y k y for fermions in regions 1 and 
2 in Fig. 1, and by e*, = —v x k x for fermions in regions 3 
and 4 in Fig. 1. 

Magnetically-mediated interaction r a 7 i ^ 5 (q, H m ) is 
proportional to dynamical spin susceptibility y(k, f2 m ): 

ka7,/3fS(k, D m ) — * @^5 X(^ (-0 

The static susceptibility y(k — K,0) comes from non- 
critical high-energy fermions (i.e., fermions with ener¬ 
gies of order bandwidth), and its momentum dependence 
is fully analytic and can approximated by a conven¬ 
tional Ornstein-Zernike form l/(k 2 + £ -2 ). The dynam¬ 
ical part (the Landau damping) comes from low-energy 
fermions and must be computed together with fermionic 
self-energy E(fc,cu m ). To simplify presentation, we fol¬ 
low earlier work s 32 : 36 ' 39 and neglect the momentum de¬ 
pendence of the self-energy, i.e., approximate E(/c,u; m ) 


by E(w m ). We comment on this approximation below. 
For E(fc,w) « E(w), self-consistent evaluation of the po¬ 
larization operator and the fermionic self-energy then 
yields&22 

x(k-K,n) SX (t.e) = t2 + 7| ‘ | + { _ 2 . (2) 

with the Landau damping coefficient 7 = Ag/^nv^) and 


= Sgn(w m )y / Wo V'|w m | + w s f - y/uJTf 


( 3 ) 


where kh is a momentum of a hot fermion, ujo = 
9 g/( 167 r), and w sf = £ -2 /7 = (7r/4)(u F ^ _ 1 ) 2 /0- This 
self-energy interpolates between Fermi liquid form at the 
smallest frequencies and quantum-critical, non-Fermi liq¬ 
uid form at larger frequencies: 

| — ^rn, for |w m | <C W s f ^ 

|sgn(w m ) 7 |u>m|u;o, for \u> m \ » w sf ’ 


where we have defined a dimensionless parameter 
A = 

47TUF 


( 5 ) 


The Green’s functions of hot fermions are [E(w m ) = 
^(^m) T ^m] ■ 


Gi(*,w) = — 


zE(w m ) - v F k , 


-,G 2 (k,u;) =-, 


zE(w m ) + v F k y ’ 

( 6 ) 

for fermions in regions 1 and 2 in Fig. 1, and 


G 3 (k,u) = 


zE(a>„ 


v F h 


-,G 4 (fc,w) = T 




v F k x 


( 7 ) 


for fermions in regions 3 and 4 in Fig. 1. 

We now comment on the approximation E(fc,w m ) « 
E(w m ). First, despite that E(fc/ I ,w m ) is parametrically 
larger than E (k—kh, 0) at A > 1 (the difference is a power 
of A, see Ref. |32), E(/c — kh, 0) is not small compared to 
v F \k — kh |, i.e., the renormalization of the Fermi veloc¬ 
ity of a hot fermion is not weak. At the same time, the 
renormalization of the Fermi velocity does not generate 
a distinction between SC and CDW-x channels and from 
this perspective is irrelevant for our consideration. Sec¬ 
ond, the coupling A does actually depend on the location 
of k^ along the Fermi surface and diverges at £ = oo only 
at a hot spot & 32 ' 37 . This momentum dependence does af¬ 
fect the values of T sc (Refs. [37i and E 2 ) and of T cdw but 
in non-crucial way, i.e., it affects the numbers but does 
not impose qualitative changes. 


III. SC AND CDW-X INSTABILITIES 

To compare SC and CDW-x instabilities, we consider 
linearized equations for SC and CDW-x order parameters 
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FIG. 2. The superconducting vertex 4> ac [Panel (a)] and the 
CDW-x vertex 4> c< j w [Panels (b) and (c)[ in the ladder approx¬ 
imation. Panel (b) shows the relation of CDW condensate 
formed by fermions at hot spots 1 and 2 with the one formed 
by fermions at hot spots 3 and 4. In turn, CDW condensate 
formed by fermions at hot spots 3 and 4 is related in a simi¬ 
lar way to that at hot spots 1 and 2, which we do not show 
here. In Panel (c) we integrated out fermions near hot spots 3 
and 4 and obtained the self-consistent equation for the CDW 
condensate formed by fermions at hot spots 1 and 2. 

and compare at what temperature these equations have 
non-trivial solutions. The equation for SC order param¬ 
eter involves conventional ladder series in the particle- 
particle channel [see Fig. 2(a)]. Each cross-section con¬ 
tains two fermionic Green’s functions with equal e*, and 
opposite frequencies and one spin-fluctuation propaga¬ 
tor. The equation for CDW-x order parameter is a 2 x 2 
set of coupled ladder equations for order parameters in 


hot regions 1-2 and 3-4. Each cross-section still contains 
the product of two G’s and one x, however in region 1- 
2 the two fermions have equal frequencies and opposite 
ek, while in region 3-4 they have equal frequencies and 
equal efc. The product of the two G’s in the region 1-2 is 
the same as the product of the two G’s in the SC cross- 
section, up to overall minus sign. The sign change is 
compensated by the summation over spin indices within 
the cross-section^: in the spin-singlet SC channel, the 
spin structure of the particle-particle vertex is ia v , and 
icr^cTap ■ a 7 s = (—3)1(7^, while in CDW-x channel the 
spin structure of the particle hole-vertex is a <5—function, 
and 5 ai <j a p • d>y = (+3)<5g,5. As a result, the kernel in 
region 1-2 in the CDW-x channel is equivalent to that in 
the SC channel. Given this equivalence, it is convenient 
to re-arrange the ladder series in the 2x2 set for CDW-x 
order parameter and represent then as a single ladder 
series involving fermions in the region 1-2 with the effec¬ 
tive interaction x'eff coming from second-order composite 
process in which fermions from the region 1-2 scatter into 
the region 3-4 and then scatter back into the region 1-2. 
This Xeff is the convolution of the two G’s in the region 3- 
4 and two bosonic \ with internal momenta/frequencies 
in the region 3-4 and external momenta/frequencies in 
the region 1-2, see Fig. 2(c). This composite Xeff has to 
be compared with the original \ a t relevant frequencies 
and momenta. 

It is instructive to consider separately the case of mod¬ 
erate £, when the pairing involves fermions with energies 
below w s f, and the case of large enough £, when relevant 
fermionic energies exceed w s f. 


A. SC and CDW-x instabilities in the Fermi-liquid 
regime 

At energies below w s f, fermionic self-energy has a Fermi 
liquid form E(w) = Xuj and the momentum and frequency 
integral of the product of the two fermionic G’s in either 
SC channel or in CDW-x channel in the region 1-2 yields 
a conventional Cooper logarithm log(w s f/T). The loga¬ 
rithm comes from the smallest fermionic u and e*,, hence, 
to logarithmic accuracy, to and ek can be set to zero in the 
interactions. The linearized equations for superconduct¬ 
ing and CDW condensates, 4> sc (fc x ) and 4> c d w (fc£c) are (in 
both cases k x is along the FS) 


J 


$ sc{k x ) = ypy log (u si /T)± f dk' x <I> sc (k' x )x(k x - k' x , 0) 

$cdw (k x ) = yipy log(w s f/T)^- f dk' x <& c a w {k' x )xcom.{k x , k' x: 0), (8) 


where A is given by ([SJ) and the shift by K is absorbed 
into the definition of x in Ed- 0- To get the overall 
sign in the r.h.s. of the equation for <f> sc (k x ) positive, 


we additionally assumed that SC order parameter has 
d-wave symmetry, in which case <I> sc (k) = — < l> S c(k + K). 
[For CDW-x channel, $ c dw(fc) has opposite sign in the re- 
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gions 1-2 and 3-4, but not equal magnitude. This implies 
that the form-factor for CDW-x order is an admixture of 
s —wave and d-wave components (an admixture of a true 
CDW order and a bond charge order), and that d-wave 
component is largei— ~— .] 

The composite Xcom(k x , k' x , 0) is the convolution of two 
dynamical spin susceptibilities and two Green functions 
of fermions with parallel velocities. Because we already 
have logWsf/T in the prefactor in (|5|). we can evaluate 
Xcom{k x , k x ,0) at T = 0, by replacing the summation 
over Matsubara frequencies by integration. 


v (i. v ni - f dn m dp x dp v 

Xcom(Kx, K x , U) - J (inm(1 + A) _ vppx) 2 

*X{K - Px,Py, tyxiK - Px,Py, ty 


where momenta p x , p y and frequency f l m are for fermions 
in region 3-4. 

Because the integrand contains double pole, a regular¬ 
ization is required. It is provided by either taking the 
external frequency for the vertex <f> c dw to be infinitesi¬ 
mally small but non-zero, or by shifting by infinitesimal 
amount the momentum (0, Q) from the distance between 
hot points. Because we analyze the emergence of the 
static CDW order and the anti-nesting between regions 
3 and 4 is only approximate , the correct regularization 
procedure is to shift the momentum. One can easily make 
sure that this is equivalent to keeping the integrand as 
in ([9]) and integrating first over p x and then over D. Be¬ 
cause the integrand vanishes at larger \p x \, the integral 
over p x over real axis can be extended in a standard way 
onto a complex plane of p x and the integration contour 
can be closed in the half-plane where there is no double 
pole. The spin-fluctuation propagator x depends on p x 
and has poles in both half-plane of complex p x . Taking 
the contributions from the poles in the two %’s in the 
half-plane where there is no double pole, and integrating 
then over p y and (in any order), we obtain^ 


Xcom(k x ,k x ,0) = x{k x - k' x ,0)A(k x £,k' x £, H^), (10) 

The function A{x, y , z ) is the scaling function of all three 
arguments and is 0 ( 1 ) when the arguments are of order 
one. When x = k x £ and y = k x £ are non-zero, A(x, y, z) 
evolves but remains close to A( 0 , 0 , z) for relevant x,y = 
0(1). As a result, to good accuracy, Xcom(k x , k' x , 0) and 
x{k x — k' x , 0 ) differ just by a constant A( 0 , 0 , z) = A(z). 
The evaluation of A(z) yields A(l) = O.llg/fa'yvp) = 
0.084, A(z 1) ~ l/(2z). At weak coupling [small A 
and hence large z = (1 + A)/A], A(z) is small, but at 
A > 1 (hence smaller z), A(z) becomes of order one. We 
plot A(^) in Fig. 3. 

Note in passing that the value of A(z) can be fur¬ 
ther increased if we abandon the self-consistent approach, 
in which 'yv'p = (4/-7T )g and assume that CDW-x or¬ 
der emerges from some pre-existing pseudogap state 
which additionally reduces the Landau damping coeffi¬ 
cien t 17 ! 44 ' 45 due to reduction of a low-energy fermionic 


Mz) 

4(1) 



FIG. 3. The plot of A(z)/A( 1) = A(0, 0, z)/A(0, 0,1) as a 
function of 2 = (1 + A)/A. 


spectral weight in the hot regions. Because A(z) is 
inversely proportional to 7 , it increases when 7 gets 
smaller. This in turn increases y C om compared with y. 

To obtain T sc and T c d w , one needs to solve Eqs. ©• 
Because 4 > sc (fc) has a finite value at k = 0 and relevant 
momentum deviations from hot spots are small (of order 
£ -1 ), one can safely approximate $ sc (fc) by $ sc (fc/ l ) and 
explicitly integrate over k x — k' x in the spin susceptibility. 
This leads to a familiar Fermi-liquid result 


1 = 


A 

1 + A 


log 


^sf 

T~ 

-*• sc. 


( 11 ) 


where u; s f differs from cj s f by a number. Hence 

l + A 

T sc — ^ ^ 


( 12 ) 


For CDW-x channel, the evaluation of T c d w requires one 
to solve the actual integral equation in momentum be¬ 
cause Xcom depends separately on k x and k x . But, like 
we said, for relevant k x £ ~ k x £ ~ 1 , A(x, y , z) can be well 
approximated by A(0, 0,z) = A{z). Using this approx¬ 
imation, we immediately find that $> c dw(k) can be also 
replaced by its value at £ 7 , and T c d w is determined from 


1 


= A- 


1 + A 


log 


^sf 

ledw 


(13) 


Hence 


i+A 

Tcdw = w s f e 


(14) 


At small A, T c d w is exponentially suppressed compared 
to T sc , but at A > 1, T c dw is only numerically but not 
parametrically smaller than T sc . 


B. SC and CDW-x instabilities at larger £ 

The Fermi liquid consideration is useful for the un¬ 
derstanding why T c dw becomes comparable to T sc at 
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A > 1, but it cannot be extended to larger £ and hence 
larger A = ?>gt; / (Airy F ) because for these £ the pair¬ 
ing comes from energies larger than w s f. Specifically, 
there are two characteristic scales in the problem: w s f, 
which is the upper boundary for Fermi-liquid behav¬ 
ior, and loo = 9<?/(167r), which is the upper bound¬ 
ary for non-Fermi liquid, quantum-critical behavior with 
E(w m ) ss sgn(w m )(|w m |a;o) 1 ^ 2 - At frequencies above loq, 
the self-energy preserves its non-Fermi liquid form but 
gets smaller than the bare u> m . The ratio of the two en¬ 
ergies is w s f/w 0 = 1/(4A 2 ) (Ref. HU). At A < 1, they 
are comparable, but at large A, w s f <C too- In this re¬ 
spect, the upper boundary of the Fermi liquid regime is 
parametrically smaller at A > 1 than the highest up 
to which S(uj m ) is relevant. In the Fermi liquid descrip¬ 
tion, the gap equation involves only frequencies ui m < lo s { 
and hence both T sc and T cc j w scale with w s f and vanish 
at A = £ = oo. Meanwhile, the quantum-critical be¬ 
havior of the system extends to wo, which remains finite 
at £ = oo. Earlier studies of d-wave superconductivity 
foun d 36 : 37 that T sc (A = oo) is in fact finite and is of or¬ 
der ujq■ The issue we consider below is whether T c dw also 
remains of order u >o at infinite £. We argue that it does 
and the ratio T c d w /T sc is larger than in the Fermi liquid 
regime. 

I 


Solving the linearized equations for & S c(k x , w m ) and 
$cdw(fci,w m ) in the quantum-critical regime is rather 
involved procedure as these equations become inte gral 
equations in_frequency for <h sc (fc x , w m ) (Refs. 3, UEIH, 
[dll. [~12l. andl46h and in both frequency and momentum for 
$cdw (k x ,u) m ). We refrain from presenting the details of 
the solution of the integral equations, but rather focus on 
proving that (i) the mean-field T cc j w is non-zero, no mat¬ 
ter what the magnitude of the composite interaction is, 
and has power-law rather than exponential dependence 
on the interaction strength, and (ii) the effective coupling 
in the CDW channel is smaller numerically but not para¬ 
metrically than that in the SC channel, hence T c and 
T c dw differ by coupling-independent numerical factor. 

To show that T c d w is non-zero for any interaction 
strength, we set T = 0 and consider the equation 
for CDW-x order parameter as eigenvalue/eigenfunction 
equation. We show that the eigenvalue E c d w is infinite at 
T = 0 for any value of the coupling. Because the transi¬ 
tion occurs when E = 1 and E decreases as temperature 
increases, the very fact that E c d w is infinite at T = 0 
implies that the instability temperature T c d w is finite. 

We first briefly demonstrate how this works for su¬ 
perconductivity. The eigenvalue equation for SC order 
parameter is 


Tsetse {k x , Cd m ) — 


35 


87T 2 Uf 


dioL 


l^ml A 


dk x ^ sc (k x ,uj rn )x(k x k x , uj m w m ) 


(15) 


To obtain this equation we used the non-Fermi-liquid 
form of the self-energy, and integrated over the transverse 
momentum k y . Because the interaction \ only depends 
on momentum transfer k x — k' x , & ac (k x , w m ) = <f , sc (u; m ) 
is a solution. Integrating over k x in the r.h.s. of Eq. m 
we obtain^ 

T sc $ sc (w m ) = - / dUJ ' m - 1 

4 J V\ U 'm\ \/|Wm - U>’ m \ 1 + yJ\u>' m \/UQ 

(16) 


One can easily verify that <f> sc (oi; m ) = 1/^/\u> m \ is an 
eigenfunction, and the corresponding eigenvalue is infi¬ 

I 


nite: 


E r , 


1 r° <**'„ _ 

2 7o Kl 


(17) 


The logarithmical divergence of E sc at T = 0 resembles 
that in the s tanda rd BCS theory, but the eigenfunction 
fkcdw ~ 1/vl^mJ is different. The divergence indicates 
that at T = 0 the normal state is unstable towards form¬ 
ing a SC condensate, and there exists a finite T sc ~ ui 0 
at which SC transition occurs. 

We now apply the same logic to the analysis of -Ecdw 
for CDW-x order parameter. The ladder equation for 
< f’cdw(w m ) in the quantum-critical regime is 


^cdw^cdw^n^m) — 


3ff 


dto' 


1 


8tt 2 df J \u>' m \ + i/l w ml w o - 

where the composite interaction in the quantum-critical regime is given by 

dfl m dp x dpy 


dk X A> c( \-w{k x ■; ^ m)x com (k x ,k x ? 1 Wm ) 5 


(18) 


3 g 


Xcom(k x , k x , UJ m , UJ m ) — „ O 


8tt 3 J (isgn(D m )(|D m |wo) 1/2 - v F p x ) 2 

1 


1 


{Px - k x ) 2 +pl + 7 I- uim I {p x - k ' x ) 2 +pl + 7|^m - <*4 


( 19 ) 
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We assume and then verify that the eigenfunction of Eq. 
m takes the form 

<f = —= 

| 


where k = k x /\J r y|w m |. In principle, ip(x,y) has both 
even and odd components in both variables. However, 
substituting ip into (fT51) we find after simple algebra that 
divergent (f dw^/lo/J) contribution to the r.h.s. of (fT51) 
comes solely from the even component <^(|fc|). In explicit 
form we have 


ip(k, sgnu m ) (20) 


1 C^° rh / f 

£ cdw ^(|fc|) = 2- / / dk'K(k,k'M\k'\), (21) 

J —ujq l^ral J 


where 

r,x __3_ [ dx dy dz ( x 2 - 9|zj/64) 

’ 327T 2 J (x 2 + 9|z|/64) 2 

1 1 

(x — k ) 2 + y 2 + |2| (x — k 1 ) 2 + y 2 + \z — 1| 

( 22 ) 

This is integral equation in k with non-singular momen¬ 
tum dependence in y(\k\). We verified that for relevant 
k < 1, </5(|fc|) can be reasonably well approximated by 
a constant. Specifically, if we substitute <f(k') = into 
the r.h.s. of (1211) we find that f(k) = f dk'K(k', k) is a 
slowly varying function of k (see Fig. 4). Taking /(0) for 
an estimate, we obtain 


E, 


cdw 


o r° dw' m 

2 Jo Kl‘ 


(23) 


where C = 0.96/7T = 0.31. We see that T cdw diverges 
logarithmically at T = 0, like E sc . As the consequence, 
at T = 0 the system is unstable towards forming a CDW- 
x condensate, hence T cdw must be finite. We note, to 
avoid misunderstanding, that at this (mean-field) level 
we consider SC and CDW-x instabilities as independent 
on each other. 

Comparing Eq. (E51) and Eq. im we see that in the 
quantum-critical regime the effective dimensionless cou¬ 
pling in the CDW-x channel is weaker than that in the 
SC channel only by the numerical factor C. We now 
use the result of the generic analysis of the SC quantum- 
critical problem^!, which shows that the dimensionless 
coupling f3 appears in the formula for the critical tem¬ 
perature T sc , as /3 2 in the overall factor, rather than in 
the exponent. Applying this also to Tcdw, we find that 
in the quantum-critical regime T cdw is smaller than T sc 
roughly by C 2 ~ 0.1. The ratio T cdw /T sc can again be 
enhanced if we assume that CDW-x emerges from a “pre¬ 
emptive” state in which Landau damping is additionally 
reduce d 17 ’ 44 ! 45 . The ratio T cdw /T sc also get enhanced 
when we include into the analysis pair-breaking effects 
by thermal fluctuations^. 


/ 



FIG. 4. the plot of f(k) as a function of k. For relevant k < 1 
it is a slowly-varying function. 


IV. THE ROLE OF A FINITE DEVIATION 
FROM NESTING/ANTINSTING AT HOT SPOTS 

We have demonstrated that for a nested/antinested 
FS, T c dw is finite for any magnetic correlation length 
only numerically lower than T sc . For a generic FS, there 
is some small but finite angle between Fermi velocities 
at, say, hot spots 1 and 2. 

We note that our analysis for the quantum-critical 
regime (large £) did not require any particular nesting 
condition on the FS - indeed, at £ = oo, for a generic 
FS one can use the same scaling form in Eq. (1201) and 
by the same scaling argument factor out f dw m /|w m | di¬ 
vergence, only the form of K(k,k') is now more compli¬ 
cated. Therefore for a generic FS the CDW-x instability 
still occurs in the quantum-critical regime, although the 
onset temperature T cdw is smaller than in the perfect 
nesting/antinesting case. At a finite £, the divergence of 
T c dw at T = 0 is, however, cut, and for any finite devi¬ 
ation from nesting/antinesting limit there exists a finite 
critical £ cr at which T cdw vanishes. 

At small deviations from nesting/antinesting critical 
£ cr and critical X cr are both small, and the computa¬ 
tion of can be done in the Fermi liquid regime. At 
small angle a = v x /v y between Fermi velocities at hot 
spots 1 and 2, the fermionic dispersion takes the form 
ei,2(k) = VF{±k y + ak x ) The convolution of the two 
Green’s functions at hot spots 1 and 2 now gives 



(24) 


where k x , y = VFk X:V /(l + A). Taking a = 0 leads us back 
to the log(w s f/T) in Eq. (7). For a/0, the logarithm 
is cut by k x . As a consequence, there exists a critical 
f — at which T c dw vanishes. Fiom Eqs. © and COD 
the typical value of k x relevant for CDW-x is £ 1 . Then 

















t; cr is given by 


V. SUMMARY 


1 = A 


X C r 

1 + X cr 


log 


fa 1 + X cr \ 
\ an X cr J 


(25) 


where a = 0(1) and we recall that A = A(z ), z = 
(1 + A)/A. Using A(z 1) « 1/(2 z), we obtain, to 
logarithmic accuracy, 


Xcr 


21oga| 1 / 2 3 4 5 6 ’ 


_ 47 TV F 

Scr 0 - 

3 9 


(26) 


Finally, we comment on the validity of the expansion 
around hot spots, which has been adopted throughout 
this work. From Eq. (9) we see that the relevant mo¬ 
menta k x ,k' x (which, we remind, are deviation from hot 
spots 1, 2 along the FS) are of order £ _1 . The typical mo¬ 
menta p x and p y in Eq. (8) are also of order £ . Taking, 

e.g., £ = 3a, we obtain that typical momentum deviation 
from a hot spot is ~ 0.06 x 2-rr/a. This is a fairly small 
momentum range. As a comparison, for the dispersion 
taken in e.g., Ref. 0, the separation between neighbor¬ 
ing hot spots is a few times higher: 0.2 x 27r/a. In this 
sense already for £ = 3a, SC and CDW instabilities come 
from the vicinity of hot spots. The approximation gets 
even better when £ increases. 


To summarize, we have shown explicitly that anti¬ 
nesting for a half of hot spots does not prevent the in¬ 
stability towards CDW-x order [the one with momentum 
(i Q , 0) or (0, Q)\ as in the strong coupling regime the cor¬ 
responding T c dw differs from T sc for d-wave superconduc¬ 
tivity only by a constant. For the case when the Fermi 
surface is not perfectly nested/antinested at hot spots, 
we found that at large £ CDW-x instability still emerges, 
but terminates at some critical £ cr ~ vf/9 - 

The ratio T c d w /T sc is still small numerically it is 
C 2 ~ 0.1 at large £ in fully self-consistent theory From 
this perspective, it is likely that spin-fluctuation ex¬ 
change is not enough and an additional mechanism, i.e., 
electron-phonon interactio n 13 ! 47 , additional softening of 
fermionic damping due to pseudogap physics separate 
from charge orde r 17 ! 44 ] 45 , or Coulomb repulsion between 
nearest neighbors^ is needed to make CDW-x a strong 
competitor to d-wave superconductivity in the pseudogap 
phase of the cuprates. Still, from theory perspective, it 
is essential that the presence of anti-nesting parts on the 
Fermi surface is not an obstacle for CDW-x instability. 
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